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Abstract. We study the simultaneous semi-classical and adiabatic asymp- 
totics for a class of (weakly) nonlinear Schrodinger equations with a fast peri- 
odic potential and a slowly varying confinement potential. A rigorous two-scale 
WKB-analysis, locally in time, is performed. The main nonlinear phenomenon 
is a modification of the Berry phase. 



1. Introduction and scaling 

In this work we study the asymptotic behavior as e — > of the following semilinear 
initial value problem (IVP): 

f ied t r = - yAV/ + V T (jj) ^ + U{xW + e\(t) m 2a i> £ , 

\r\ t=0 = i> e Ax), 

where a; e l* 1 , t e I, c e N and < £ < 1. Here and in the following e- 
dependence will be denoted by the superscript e. The external (confining) potential 
U = U(x) £ K is assumed to be smooth on R d , whereas the lattice-potential 
Vt = Vr(y) £ R is assumed to be smooth, uniformly bounded in R d and periodic 
with respect to some regular lattice T ~ Z d , generated through a basis {£i, . . . , Cd}, 
Q £ R d , i.e. 

(1.2) Vr(y + ~f) = V r (y), V2/eR d , 7 er, 
where 

(1.3) r=| 7 eR d : 7 = ^7/0, 7/ezj. 

Finally, we assume A = X(t) £ K to be a smooth coupling-function and x/jf £ L 2 (R d ) 
to be normalized such that 



(1.4) / \xl) £ I {x)\ z dx = l. 

This normalization is henceforth preserved by the evolution since X(t) £ R. 

Nonlinear Schrodinger equations (NLS) of type 1|1.1|) appear in various physical 
situations, cf. [44] for a general overview. An important example in d = 3 is the 
case (7 = 1, X(t) = ±1, i.e. the so called repulsive resp. attractive Gross-Pitaevskii 
equation, a celebrated model for the description of the evolution of Bose-Einstein 
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condensates (BECs) 36 . In order to motivate the scaling in Ijl.lJ) we shall examine 
this case more closely: 

In physical units, the Gross-Pitaevskii equation (for d = 3) is given by |36| 

h 2 

(1.5) ihd t ^ = — - — + V fx) ip + U (xU ± Na(t)\ip\ 2 ip, 

2m 

where m is the atomic mass, H is the Planck constant, N is the number of atoms 
in the condensate and 

(1.6) a(t) = , 

with a(t) € R denoting the s-wave scattering length derived from the corresponding 
iV-particle theory, cf. jSOlEEl- (The fact that a(t) is chosen time-dependent is 
motivated by recent experiments on BEC where this has indeed be achieved by 
some highly sophisticated experimental techniques.) In this context the external 
potential U(x), which traps the condensate, is usually assumed to be a harmonic 
confinement potential of the following form [2 ^] : 



mui 



(1.7) U (x) = ^\x\ 2 , lu eR,xeR*. 

More general, non-isotropic variants of such confinement potentials are used to cre- 
ate so called disc-shaped or cigar-shaped, i.e. quasi two or, resp., one dimensional, 
BECs (see |2 EU and the references given therein) . If in addition a periodic po- 
tential V(x), which in physical experiments is generated by an intense laser field, 
is included, the condensates are referred to as lattice BECs. A particular example 
of V is then given by 

(1.8) V{x) = ^^-^?^m), 

i=i 

where £ = (£i, £2, £3) with ^ £l denotes the wave vector of the laser field The 
sign in front of the nonlinearity in 11.5fl corresponds to a stable (defocusing) resp. 
unstable (focusing) condensate. To rewrite the equation (|1.5(l into our semi-classical 
scaling we proceed similar to 2 . More precisely, we introduce dimensionless vari- 
ables 

(1.9) t = w t, x = —, i>(t, x) = x 3 / 2 i)(t, x), 

X s 

where x s will be determined later and ip(t, x) is such that the normalization (|1.4(l 
is preserved for d = 3. Multiplying 1)1. 5fl by l/(mw 2 x 2 ) and omitting again all " ~" 
we find the following dimensionless equation: 

2 

(1.10) iedti) = -y A ^ + ( j) V> + U(x)tp ± S(t)e 5/2 \ip\ 2 ^, 
where the potentials are defined by 

V{xsey) 

and the appearing parameters e, S(t) € R+ are 



(1.11) V T (y) := v | f , U{x) := ^- 

mujAxi 2 



(1.12) e :=-A_ = f5aV, 6(t) := ^ = 47r|a( * )|iV 

toomxj \x s J a^nuo ao 
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with ao denoting the length of the harmonic oscillator ground state corresponding 
to Uq(x), i.e. 

(1.13) ao :=J—. 

V uJom 

Since we aim for e -C 1 and Se 5 ^ 2 to be of the order of e we require 8 = 0(e~ 3 / 2 ), 
hence 47r|a|iV ^> ao, which from a physical point of view corresponds to the strong 
interaction regime, also known as Thomas-Fermi regime )36j . Now, consider a 
reference value a for a(t) and similarly denote by 6 the parameter 6 for this reference 
value a. Inserting 1)1.12)1 into Se 5 ^ 2 — e, we compute the characteristic length scale 



(1.14) x s = (AnNldla 2 ) 1 ^, 

which one needs to choose as the appropriate reference scale in our situation. In 
particular we shall assume \ipj (x)\ to vary on this scale. The coupling function X(t) 
is then given by A(i) = S(t)/S. Identity (j 1 . 1 4|l implies 

< L15 > -(i^kf* 1 ' 

which is different from the one given in Moreover, having in mind 1)1.8)) . (|1.11|) 
we require for the periodic potential Vr 

(1.16) e£ix s = 0(l), \% 2 = 0(1), for/ = 1,2,3. 

ZiTfl X s UJq 

From these relations one computes 

(1.17) & « a - 4/3 (4^|a|) 1 /3, for /= 1,2,3, 

which gives the required wave vector in our regime and one checks that in this case 
the conditions 1)1.16)1 are satisfied. We remark that this scaling is in good agreement 
with some typical recent experiments. For example in the case of a lattice BEC 
consisting of Rb atoms we have, cf. |21 111): 

(1.18) a «3,4xlO" 6 H, a « 5, 4 x 10~ 9 [toL A^l^xlO 5 . 

This gives: 47r|a|iV « 10 _2 [m] 3> ao, hence e»4.3x 10~ 3 -C 1 and for the wave 
vectors we compute w 4, 6 x 10 6 [l/m], which is of the same order of magnitude 
as stated in [§]. The reference length scale in this case is x s — 2. 1 x 10 _6 [m], which 
is O(ao). Finally, to motivate the choice a > 1, we note that for d < 3 higher order 
nonlincarities are frequently used in the description of BECs [351 HID] • 

From a mathematical point of view the limit e — > corresponds to the simultaneous 
semi-classical (or high-frequency) and adiabatic limit (see [32, 40, 45 for general in- 
troductions to these fields). For linear time-dependent Schrodinger equations (with 
periodic potentials) this asymptotic regime has been intensively studied by several 
authors, using (spatial) adiabatic decoupling theory |341 145) or Wigner measures 
J! El UHl- to mention results obtained in recent years. A numerical study of these 
asymptotics can be found in 21 . 

In our scaling the nonlinearity is o(l) and can thus be called weak, still it makes the 
rigorous asymptotic analysis of the given I VP considerably harder. Even without a 
periodic potential the semi-classical limit for NLS is still far from being completely 
understood. In particular, we cannot use the above mentioned mathematical tech- 
niques, which so far only work in a linear setting. (For a notable exception see 0.) 
Thus we shall rather apply a more naive asymptotic expansion method in the spirit 
of the traditional WKB-type expansions. Due to the periodic potential, we use a 
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so called two-scale WKB-ansatz, first introduced in 4 , which has already been 
successfully applied in the case of linear periodic Schrodinger equations ^1 EHj • 
Our scaling is such that the nonlinearity enters in the leading order term of the 
asymptotic WKB-type solutions, although the Hamilton- Jacobi equation for the 
phase of the wave-function is found to be the same as in the linear case. This is 
analogous to the weakly nonlinear ( dispersive ) geometrical optics regime discussed 
in . (See also 02] for an application of this scaling in another semi-classical 
context). The asymptotic description is valid on macroscopic time-scales t = 0(1) 
but in general only for small \t\ > 0. 

Before giving a precise description, we state the typical result that we shall prove. 
The possibly not well-defined assumptions in the following statement will be dis- 
cussed more precisely below. 

Theorem 1.1. Let d > 1, and U be smooth, real-valued potentials, Vp being 
r -periodic, U being sub- quadratic, and A being real-valued and smooth. Assume 
that the initial datum ipj is of the form 

riix) = ai (x) X n (p V&(a;)) e*+*W e +e<fi(x), 

where ai £ 5(R d ;C), <f>i G C°°(R d ;R) and \n = Xn(y,k) is a Bloch eigenfunction 
associated to a simple isolated Bloch band E n = E n (k). We suppose that ip\ satisfies 
Assumr>tion \3.5\ below, with K > d. Assume that no caustic is formed before time 
t > 0, and fix r G]0, t[. Then there exists £o > such that for < e < Eq, the 
solution t[) £ to is defined up to time r . Moreover, it satisfies the following 

asymptotics as e — > 0: 

sup W e (t)-vl(t)\\ L2m =0(e), 

0<t<T O 

(1.19) sup H e (t)-Vo(t)\\L°om«) ^0(e 1 -^), for any r, > , 

0<4<T 

where the approximate solution Vq is given by: 

vS(t,*) = Xn (*,V x <t>{t,x)) eHVTWeW'. 

Here, <fi solves the Hamilton- Jacobi equation (|2.9I) . corresponding to the classical 
flow: (t,x) i— » X t (x), as defined by (|2.15|l . Jt is the associated Jacobi determinant 
(I2.1ti|l , and to is given by 

(j(t,x)=—i / (3 (s, X s (x)) ds 
Jo 

-M*)| 2ff f I \Xn(y,V x 4( S ,X s (x)))\ 2 ° +2 dyds. 

Jo ■Js(x) Jy 

We denote by (3 6 iM. the Berry phase (|3.6(l . and by Y the centered fundamental 
domain of V . 

Remark 1.2. Our result holds only before caustics. This should not be surprising; 
even in the linear case A = 0, the WKB method is effective only away from caustics. 
On the other hand, some techniques have proved to be efficient to overcome this 
difficulty in a linear framework, such as Gaussian beams (see e.g. ^2]) or Wigner 
functions (see e.g. [21 ESI)- However, adapting these techniques to a nonlinear 
context seems to be a challenging open question. 
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Remark 1.3. The assumptions on the corrector ipj for the initial data are not 
trivial (see Assumption 13 . 5|1 . They state essentially that the initial data are well- 
prepared, in order to prove a nonlinear stability result. Note however that ifj is 
of order O(l)ase^0in any reasonable sense. The assumptions K > d means 
that we have to consider (at least) d correctors to prepare the initial data. This 
assumption may seem surprising; the proofs we give rely on it, and it would be 
interesting to understand how necessary this assumption is. 

The above result shows that the leading order nonlinear phenomenon is represented 
by the phase factor uj. The Berry phase is a linear (geometrical) feature (see ()3.6(l 
below), but the second integral in the definition of uj stems from the nonlinearity. 
In the context of laser physics, this phenomenon is known as phase self- modulation 
(see e.g. E3E1[T3]). 

The paper is organized as follows. In Section we start a formal asymptotic 
expansion, following WKB-methods. This leads us to consider the Bloch eigenvalue 
problem. The asymptotic expansion is considered in more detail in Section|3| where 
a formal approximate solution is constructed at any order. The justification of this 
approximation is performed in Section We discuss our results and some of their 
possible extensions in Sectional In Appendix lAl we detail a computational step 
from Section 

2. Asymptotic expansion: emergence of Bloch bands 
For solutions of (|l.lfl we seek an asymptotic expansion of the following form: 

oo 

(2.1) ^(t,x)= u £ (t,x,-)e l ^' £ ; u e {t,x,y)~ £Vu 3 -(t ja; ,i/), 

where we assume that both (j)(t,x) S R and u E (t,x,y) S C are sufficiently smooth. 
Moreover we impose 

u e (;-,y + j) = u e (;-,y) 1 V2,eK d , 7 er. 

We assume that the initial condition is compatible with (|2.1|) : 

Assumption 2.1. The initial wave-function ipf is in the Schwartz space S(M. d ), 
and is of WKB-type, i.e. 

(2.2) ^j(x) = uj (x, ^) e tf'(*)/" + e <p°(x), 

with 4>i e C*°°(M d ;M), m £ S(R d x T d ;C)\ T d = R d /T. The function <pf is a 
corrector to be precised later on. 

From now on we shall denote the linear part of the Hamiltonian operator by 

(2.3) H * :=~/± + V T (^)+U{x) 
Plugging the ansatz (|2.1|l into Ijl.lfl we (formally) obtain: 

ied t ^ £ - H £ i]f - e\(t)\ip £ \ 2r7 i> £ = b £ (t,x, ^ e l ^ x)/e . 
We consequently expand the r.h.s. of this equation as 

oo 

(2.4) b £ {t,x,y)^Y. £lb ^ x ^v) 

j=o 

^That is, uj is rapidly decaying w.r.t. the first variable (x), smooth w.r.t. the second one (y). 
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and choose the asymptotic amplitudes Uj in a way such that bj(t, x, y) = 0, Vj > 0. 
Setting b (t, x,j) = yields 

(2.5) - iV x cj> ■ S7 y u + + V r (y)u + (U(x) + d t (j>) u Q \ y=f = 0. 

Uncorrelating the variables x and y, we shall seek a solution to the more general 
equation: 

(2.6) -^p - iV x <j> ■ V y u a + ^^-u + V r (y)u = - (U(x) + d t 4>) u . 
Denoting by 

(2.7) H r (k):=^(-iV y + k) 2 + V T {y), k £ R d , 
we can rewrite equation l|2.6[) in the following form: 

(2.8) HHV^Juo = - (U(x) + d t (j>) u . 
We now require that for some fixed n £ N, it holds 

(2.9) E n (V x cf>) = - (U(x) + 8t<t>) , 

where E n (k), k £ M. d , is the n-th eigenvalue of the Bloch eigenvalue problem |S]: 
f H r (k) X n{y,k) = E n (k) X n{y,k), n£N,y£Y, 
I Xn(y + l,k) = Xn(y,k), for 7 £ r. 



(2.10) 



Here and in the following, we denote by Y the centered fundamental domain of the 
lattice r, i.e. 



1 1 

2' 2 



(2.11) Y :=t 7 eR d : 7 = ^7i0, 7/6 

whereas Y* , denotes the corresponding basic cell of the dual lattice T* . In solid 
state physics Y* is called the Brillouin zone hence we shall denote it by B = Y*. 
Let us recall some well known facts for this eigenvalue problem, cf. |331 1451 l4l)| : 

Since Vr is smooth and periodic, we get that, for every fixed k £ B, H-p(k) is self- 
adjoint on H 2 (T d ) with compact resolvent. Hence the spectrum of Hr(k) is given 

by 

a(H r (k)) = {E n {k) ; n £ N*}, E n {k) £ R. 

In general we can order the eigenvalues E n (k) according to their magnitude and 
multiplicity, 

Ex(k) <-..<E n (k) <E n+1 (k) < ... 

Moreover every E n (k) is periodic w.r.t. T* and it holds that E n (k) = E n (—k). The 
set {E n (k); k £ B} is called the nth-energy band. The associated eigenfunction, the 
Bloch waves, Xn{y, k) form (for every fixed k £ B) a complete orthonormal basis in 
L 2 (Y) and are smooth w.r.t. y £ Y. We choose the usual normalization 



(2.12) (Xn(-,k),Xm(-,k)) L 2( Y ) = / Xn{y,k)x m {y,k)dy = 8 nt m, n, m £ N. 

Concerning the dependence on k £ B, it has been shown 125] that for any n £ N 
there exists a closed subset U C B such that: E n (k) are analytic, k) can be 
chosen to be analytic functions for all k £ f2 := B\U, and 

(2.13) £ n _i < £7„(fc) < S n+1 (fc), Vfc g fi. 
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If this condition holds for all k G B then E n (k) is called an isolated Block band |4"5] . 
Moreover, it is known that 

measW = meas{fc G 2? | E n {k) = E m (k) 7 n ^ m} = 0. 

In this set of measure zero one encounters so called band crossings. 

Equation (|2.9(l is called the n-th band Hamilton- J acobi equation corresponding to 
the semi- classical band Hamiltonian 

(2.14) h s n c (k,x) := E n {k) + U(x), (k, x) G T* x R d , 

with an effective kinetic energy given by the n-th eigenvalue for k G T* = R d /T*. 
The characteristic differential equations corresponding to (|2.9|l are consequently 
given by the equations of motion: 

, (x=V k E n (k), x\ ._= x e R rf , 

(2.15) <^ . '*- 

[ k = -V x U(x), k\ t=Q = V x (j>i(xo). 

This system (locally) defines a flow map (x, t) i— » X t (x) = X t (x; V x 4>i(x)) in physi- 
cal space. In general caustics will appear in this flow, which prohibits the existence 
of globally defined smooth solutions for Ij2.9|l . Let us denote by 

(2.16) J t {x) :=Aet{^ x X t {x;V x (j)i{x))) 

the corresponding Jacobi determinant. We have Jq(x) = 1. Denote by r the time 
at which the first caustic appears, i.e. 

(2.17) t := inf{* > | 3 x G R d : J t {x) = 0}. 

We thus have Jt(x) > for < t < r. Standard theory implies the following: 

Lemma 2.2. If h s n c (k,x) G C°°(T* x 0/ G C°°(R d ) ; i/ien there exist r > 
anrf a unique smooth solution <fi G C°°([0, r[xR d ) of the Hamilton- Jacobi equation 

M + h^iV^x) = ; 4 =Q = <Ma). 

To make sure that E n (k) (and hence hff(k,x)) is sufficiently smooth, we shall 
impose the following assumption: 

Assumption 2.3. The amplitude ui(x,y) is assumed to be concentrated in a single 
isolated Block band E n {k) corresponding to a simple eigenvalue of Hr{k), i.e. 

(2.18) ui(x, y) = a I (x)xn(y, V x (j)i{x)), 
where aj G iS(M d ;C) is a given initial amplitude. 

From (|2.8(l and (|2.10|l we conclude that there exists do = ao(t,x) such that 

(2.19) u (t, x, y) = a (t, x)x n (y, V x (f>(t, x)). 

Remark 2.4. Note that also in the linear case, assumptions similar to Assump- 
tion EI] are usually imposed, cf. |23ls2]. There however, the reason is largely to 
avoid band crossings in order to obtain global-in-time results. (The rigorous study 
of band crossings is quite involved and up to now established only for certain model 
problems, cf. [13 CHI EH • ) 

Due to caustics (and possibly additional nonlinear effects if X(t) is not real-valued, 
see Sect. EJ, we cannot hope for such global- in-time results in our case. As- 
sumption 12.31 therefore is only imposed for regularity reasons and could be sig- 
nificantly weakened, since, with some technical effort, one could modify the sub- 
sequent analysis. Indeed, all statements could be formulated locally in regions 
U C R t x M. d which neither contain caustics nor band crossings (in the sense that 
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E n {V x ^(t,x)) ^ E m {y x <j){t, x)), for all (t,x) £ U). In this way one could include 
also non-isolated bands E n (k). 

We further remark that in the case d = 1 all band crossings can be removed through 
a proper analytic continuation of the bands, cf. |39j . 



3. Derivation of the transport equations 

To characterize the principal amplitude do, we set b\ — in (|2.4|l . which yields 

(3.1) H T {V x <t>)ui + (U(x) + dt<f>) ui = L lUo - \{t)\uo\ 2a u ai 
where the linear differential operator L\ applied to uq reads 

(3.2) Liu := id t u + iV x <j> ■ V x u + i^^-UQ + div x V y u . 

We multiply equation l|3.1|l with Xn(y^x4>) and integrate over the fundamental 
domain Y. From (|2.9(l . the left hand side of 1)3.1(1 is (Hp — E n )ux; since Hr is 
self-adjoint, the integral obtained from the l.h.s. of l(3.1|l is identically zero, hence: 



(3.3) J x n (y, V x ^) (Liu - \(t)\u \ 2a u ) dy = 0, 

is a necessary condition such that (|3.1|) can be solved for u\ in terms of uq. This 
condition is known to be sufficient, from the orthogonal decomposition method 
(also known as "Feschbach method"), since E n is an isolated eigenvalue. After 
some lengthy computations, given in the appendix, we find that (|3.3() is equivalent 
to the following nonlinear transport equation for a$: 

{d t a + Ca - /3(t,x)a = in(t, x)\a \ 2c ' a , 
a o| t=0 = a i( x )- 



(3.4) 

Here, C is the usual (geometrical optics) transport operator associated to h^(k,x): 



(3.5) Ca := V ' k E n (V ' x <j>) ■ V x a + i div x (Vfc-B„(V x ^))a . 

Moreover, we have 

0(t,x) := (xn(; V x $, V kX n(; Vx<f>)) L 2 (Y) ■ S7 x U(x) 

= £<Xn(-,V B 0), ^ Xn ^ V ^L HY) -irU(x) 
i=i 

and 

(3.7) K (t, x) := -A(t) / \ Xn (y, V a 0(t, x))\ 2a+2 dy. 



This term can be interpreted as an effective coupling of the self-interaction within 
the nth-energy band. 

Note that (|2~T2*|) implies 

Re ( Xn (-, k), VfcX„(-, k)) L2[Y) = 0. 

Hence, f3(t, x) — ilm(3(t, x) only contributes a variation in the phase of do, the so 
called Berry phase 02 ESI- ^ is due to the interaction of the lattice and the slowly 
varying potential U . In our case the Berry phase in addition gets modulated in a 
nonlinear way by the right hand side of l|3.4|) . 
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Remark 3.1. The term — i/3 =: A n can be interpreted as a gauge potential, i.e. 
a connection in the (complex) eigenspace-bundle corresponding to E n (k), cf. |45j . 
For some particular lattice configurations (if the crystal has a center of inversion, or 
some other special symmetry), the curvature of the Berry connection f2„ := V x A n 
is identically zero, and the Berry connection is a closed 1-form, cf. j3SJ 1411 fH)] for 
a broader discussion on this. 

Remark 3.2. We provide a link with some already existing results. In |34l 145) the 
authors, roughly speaking, prove that in each isolated Bloch band E n (k) the linear 
Hamiltonian H £ , defined in (|2.3p . can be unitarily mapped into an effective band 
Hamiltonian h e n , which is the Weyl quantization of the semi-classical symbol 

h e n (k,x) ~ h s n c (k,x) +eh 1 (k,x) +0(e 2 ). 

This is done by constructing an e-dependent unitary operator, which block-diagonalizes 
the Bloch-Floquet Hamiltonian of the system, such that the relevant band decouples 
from the rest of the spectrum. Above the principal symbol hff(k,x) is defined as 
in (|2.14(l and the first order correction is such that 

hi(V x <f>(t,x),x) = i(3(t,x). 

Additional terms appear in h\(k, x) if one includes external magnetic fields too, cf. 

The following lemma proves that lj3.4|l has a smooth solution up to caustics: 

Lemma 3.3. Assume (f> G C°°([0,T[xR d ), and aj G 5(R d ; C). Then along the flow 
(t,x) i— > X t (x), <|3.4fl has a unique solution ao G C°°([0,T[;S(M. d )), given by: 

ai(x) ( M2(T f* k(s,X s (x)) 



ao(t,X t (x)) = exp i\ ai {x)\ 2 ° / ' \\J' ds+ / /? (s, A s (x)) ds 

Proof. Using Liouville's formula, 

— Jt{x)=&iv x (y h E n (y x <t>{t,X t {x))))j t {x) ; J (x)= 1, 

we rewrite the transport equation (|3.4|) as an ordinary differential equation along 
the flow defined by the dynamical system (|2.15|) . Let ao(t,x) := ao(t,X t ): 

' d (y/J t (x)a Q ) = /3(t,X t )a + in(t,X t ) |a | 2<T a , \t\ < r. 



y/Mxjdt 



If we define ao := y Jt (x)ao, then the principal amplitude is determined by 
(3.8) 



— a =/3 (i, <S + ik (t, X t (x)) ° d , |t| < r, 

at \Jt\x)\ 



. "°lt=o = a/(ic). 
This implies (since (3(t,x) G and n(t,x) G K) 

4- |a (t, x)| 2 = , hence |o:o(£, x)| = |a/(x)| , Vt G [0, r[. 

Define the phase shift g of <5o by ao(i, x) = ai(x)e l9 ^ t,x K Then g solves 

- g(t, x) - Im/3 (t, X t (x)) + k (t, X t (x)) 1 °^ ' ^ , 

with g\ t _ Q = 0. Inserting |oo(£) x)\ = \aj(x)\ yields the lemma, since x t—* X t (x) is 
a diffeomorphism of R d for fixed t G [0, r[. □ 
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Remark 3.4. Note that along the flow 

(t, X t (x)) = ( Xn (; V x <t> (t, X t {x))) , j t xn (; V x <f> (t, X t (x))) ) J . 3(y) , 

which is exactly the same expression as given in |23| , there however the authors do 
not distinguish between ag and ctg. 

So far we explicitly constructed an approximate solution, which solves up to 
terms of order O(e), since u\ is not fully defined yet. To obtain a better approxi- 
mation we need to set the term b 2 in l|2.4J) equal to zero, which gives 

Hr(Vx4>)u2 + (U (x) + d t 4>) u 2 = L\U\ + L 2 u - 

(3.9) 



-A(t)((2<T + l)|u | 2<T ui + 2a\u \ 2a ~ 2 u 2 u 1 



where for Uo(t, x, y) = ao(t, x)xn(y, ^x4>) we define 

L 2 u := - A^uo. 

Introduce the notations 

(3.10) L (t, x) = ffr(Vx^) + U(x) + d t cj>(t, x) 



F{z) 



\z?°z. 



From (|2.7|l . Lo is a (t, x)-dependent operator in y, and since a G N, F is smooth. 
The following projector was used to derive the transport equation (|3.4|) : 

(3.11) II„(i,ar) (^Oj(*,x)xj (y, V x (j>(t,x))^ = a n (t, x)x n (y, V x (j>(t, x)) . 

Define Q(t, x) = Id — H n (t, x). This operator is smooth, and a partial inverse for Lq 
can be defined on its range (by elliptic inversion): Lq Q is well-defined, and smooth 
(up to caustics). Applying the operator II n to (|3.9|l . the solvability condition reads 



s=0 



dy = Q. 



(3.12) J x n (|/,V x 0) (Liui+L 2 uo-\(t)-^F(u + su 1 ] 
We decompose m as 

(3.13) ux(t,x,y) = a 1 (t,x)xn (y,V x (p(t,x)) +u^(t,x,y), 
where a± is some yet unknown function and is such that 

U n (t,x)ui(t,x r ) = (Xn(;V x $),ui(t,x,-)) L 2 {Y) = 0, V(t,a:) £ [0,r[xM d . 
Now, is determined by l|3.1|l : 

(3.14) u± = L„ x g (iiuo - A(t)F(«o)) . 

which implies uj^ £ C°°([0, t[; iS(K d )), since ito is, by Lemma [3.31 Note that this 
relations imposes a particular form for the initial perturbation ipj, that is 



(3.15) Q(0, x)<fi{x) = (L^Q) (0, x) (Lwi - A(0)F(«j)) + 0(e) . 

The term 0(e) will be defined more precisely later on. On the other hand, plugging 
(|3.13|l into H3.12f> yields an inhomogcncous linear version of the transport equation 
(|3.4|) for ai (the propagating part of Ui): 



d t cii + Cai — /3(V x 4>, x)ai + i\(t)—F(uo + sai) 

as 



S(t,x), 
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where we may choose ai| t=0 = 0. The complex- valued source term S(t, x) is given 

by 

(3.16) S(i, x) = ill n {t, x) (-Liitf + L 2 u ) = i (x n {; V x 0), Liwf + L 2 u Q ) L2 . 

By this procedure, all higher order terms Uj(t,x,y), j > 1, of the asymptotic 
solution (|2.1|l can be obtained (recall that the nonlinearity F is smooth). Clearly 
we have that Uj € C°°([0, r[; iSQR. )) for all j > 1. At each step however, an 
additional condition must be imposed recursively for the initial datum ipj. This 
approach is very similar to the one followed in |14) . except that the Fourier modes 
are replaced by "Bloch modes" . 

Under the assumption (|2.ip . I|2.3|) . we construct an approximate solution, which 
solves <|1.1[) up to a remainder 0(e°°), provided that the initial data are well- 
prepared. To state precisely this property, define, for TV > 0, 



(3.17) v° N (t,x):= if N (t,x^e»W*= [±e^(t,x^))e^ 




We will use the following spaces, for s 6 N: let 

\a\ + \0\<s 

We define XI as: 



K := \f e e L 2 (R d ) ; sup ||/ E || X| < 

I 0<e<l 



These spaces are reminiscent of the spaces H^(R d ) introduced in [22] (see also 
[37]). There the dependence upon e is to recall that exactly one negative power of 
£ appears every time the approximate wave-function is differentiated. In our case, 
such negative powers also appear because of the variable y and the substitution 
y = x/e. The control of the momenta is needed because of the potential U (it 
would not be needed in the proof of Theorem 14.51 below with U sub- linear) . We 
can now state precisely the assumptions on the initial data: 

Assumption 3.5 (Well-prepared initial data). The initial data ipj satisfy Assump- 
tions \2.1\ and \2.3l and for some K £ N, the perturbation Lp\ is of the form 

K 



(3.18) <p* x (x) = e*'< a >/ s Xy-Vi(aM0 



3=1 



y=x/e 



where the O (s K ) holds in X^ for any s 6 N. The function e l ^'/ e (p\ is given by the 
first term of the right-hand side of (|3.15() . and if we denote ipa = ui, cfj(x,y) is 
given recursively for < j < K — 2 by 



= (L 1 Q) (0, x) j ii^+i + L 2 <fj - X^-^^Ffu! + 



e=i 



In the case K = 0, the sum in (|3.18(l 



is zero. 



Remark 3.6. We chose to impose II„(0, x)(fj(x, •) = for j > 1 (when we picked 
a\ | t _ Q = for instance). Our approach would also work with non-zero, smooth data 
(<Pj)i<j<K not necessarily satisfying this polarization property. All this approach 
is very similar to the one followed in |26| to justify nonlinear geometric optics for 
hyperbolic equations (see also [37], and J3] for the dispersive case). 
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We have the following Borel type lemma (see e.g. |37j): 

Lemma 3.7. There exists ipf G S(M. ) satisfying Assumvtion \3.5V such that (|3.18|) 
holds for any K G N. 

First, we will justify the asymptotics when the initial datum is given by the above 
lemma. We will then show how to relax this assumption. Note that the above 
approach is a nonlinear analog to the procedure followed in |34| . In |34| , the authors 
construct e-dependent "super-adiabatic" subspaces, in order to prove higher order 
asymptotics in the linear case. In the present context, high order asymptotics are 
needed to control the nonlinear terms (see the proof of Theorem 14.5(1 . 

Proposition 3.8. Let tpj as in Lemma \3.1\ Let r > be the time at which the 
first caustic is formed (if any). Then for any N G N, solves 



where H £ is defined by and r £ N G C°°([0,r[; S(R d )), p% G S(R d ) are such 

that r% G L™ c ([0,t[;X s £ ) and \\p e N \\x S = 0(1) for any s G N. 



To prove that the above WKB-method yields a good approximation of the exact 
solution, a nonlinear stability result is needed. First, we make our assumptions on 
the potentials precise, and establish an existence result for 1)1.1(1 . Next, we prove 
the validity of the approximation derived above. 

Assumption 4.1. The potentials are smooth, real-valued: Vr,U G C°°(IR d ;M). 

(i) Vr is V -periodic, i.e. it satisfies 1(1.2(1 . 

(ii) U is sub-quadratic: d a U G L°°(R d ) , Va G N d such that \a\ > 2. 

Remark 4.2. The assumptions on U include the cases of an isotropic harmonic 
potential (U(x) = \x\ 2 ), and of an anisotropic harmonic potential (U(x) = Yl^j^)- 
It may also be taken equal to zero, or incorporate a linear component E-x, modeling 
a constant electric field (Stark effect, see e.g. |10p. 

4.1. Existence of solutions to ((1.1|) . 

Lemma 4.3. Let Assumption ^ . l\ be satisfied, and let ipf G 5(1^), the Schwartz 
space. Lets > d/2. Then there existst e > and a uniqueip 6 G C(] — t e , t e [; H s (M. d )) 
solution to dTTTJ. Moreover, x a ip E G CQ - t 6 , f e [; H s (R d )) for any a G N d , s G N, 
and the following conservation holds: 



Proof. Since the dependence upon e is irrelevant at this stage, the above statement 
follows from the study of 



(3.19) 




4. Nonlinear stability of the approximate solution 



ill^(*)IU» = o. 



(4.1) id t iP= --^ + W{x)^ + \{t)\^\ 2 ^ ■ ip\ t=Q = ipi(x) 



where: 



• The potential W is smooth, real-valued and sub-quadratic. 

• X(t) is a smooth real-valued function. 

• a G N. 

• i>! G S(R d ). 
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The dependence of the local existence time t E upon e appears with scaling. Notice 
that the nonlinearity z i— > \z\ 2a z is smooth, because a E N. Since W is sub- 
quadratic, the Hamiltonian — + W is essentially self-adjoint on C^ (]R d ) (see 
for instance EU). The assumption s > d/2 yields H s (R d ) C i°°(IR d ). Therefore, 
local existence and uniqueness in H s (M. d ) follow from a fixed point argument, using 
Schauder's lemma (see e.g. |51ET7|'). 

To prove higher order regularity of -0 and its momenta, one can follow the proof 
of [23 (see also That article is for the case W = 0; the proof uses Strichartz 
inequalities, following from dispersion estimates. When W is smooth, real-valued 
and sub-quadratic, the same dispersion estimates are available ( |17M18p . and they 
imply the same Strichartz inequalities (|28|). Another difference with is that 
the Galilean operator x + it\7 x commutes with idt + ^A, but in general not with 
idt + ^ A — W . This is not a problem in view of the above result, since 

[x + UW X , W] = itVW = O (1 + \x\) . 

Thus, ip, xip and V x ip solve a coupled, closed system of Schrodinger equations. A 
similar argument allows to treat higher order momenta and derivatives. 

The conservation of the i 2 -norm follows from standard arguments (see 8 ). □ 

Remark 4.4. One cannot expect global existence in general. For instance, if X(t) 
is a negative constant and if a > 2/d, finite time blow-up may occur (see e.g. |5]). 
On the other hand, we shall prove below that the solution i\) z cannot blow-up before 
a caustic is formed, at least for e sufficiently small. 

Notation. Let (a e ) < £ <i and (/3 £ )o< E <i be two families of positive numbers. In 
the following we shall write 

if there exists a C > 0, independent of e e]0, 1], such that 

a £ <C(3 e , for all e e] 0,1]. 
(The C may very well depend on other parameters). 

4.2. Accuracy of the approximation. The main result we shall prove is the 
following: 

Theorem 4.5 (Stability result). Let ipj — ipj as in Lemma \3. r l\ Let r > given 
by (I2.17() . and v e N given by (|3.17() . Then for any tq s]0, t[, there exists Eq > such 
that for < e < so, the solution ip £ to (| 1 - 1 1> is defined up to time To- Moreover, for 
any JVeN and seN, 

(4.2) sup \\^(t)~v £ N (t)\\ XB =0(e N+1 ). 

0<t<T 

Proof. For N E N, we define the error term as := ip £ ~ v e N . From and 
(l3~H)j) . it solves 



(4.3) 



' ied t v% = H s v% + e\(t) {m 2 ^ e - \v e N \ 2 °v%) ~ e N+1 r% , 



e I _ N+l e 
"JVlt=0 _ ? N ' 



where H e is defined by l|2.3|l . We start with the standard energy estimate for 
Schrodinger equations: multiply the above equation by w^,, integrate over M. d and 
take the imaginary part. Since H E is self-adjoint, this yields 



edt h%(t)\\ L * < e\X(t)\ \\m 2 °r ~ k^| 2 ^|| i2 + e N+1 \\r%(t)\\ 



L 2 ■ 
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Since we work on the fixed, finite interval t G [0, to], the smooth function A is 
bounded, and the above estimate implies: 

(4.4) d t U N (t)\\ L2 < \\\^\ 2a r - W%\ 2 °v%\\ L2 +e N \\r%(t)\\ L2 . 

The idea is now to factor out in the right hand side of the above inequality, and 
take advantage of the smallness of the source term. To carry out this argument, 
we follow the method used to justify (nonlinear) geometric optics for hyperbolic 
systems; we refer to |37j for an expository presentation. 

Following 37, Lemma 8.1] we have the following Moser-type lemma: 

Lemma 4.6. Let R > 0, s G N, and F(z) = \z\ 2a z for a G N. Then there exists 
C = C(R, s, a, d) such that ifv satisfies 

\\x a (sdfv\\ Loa{Md) < R for all |a| + \f3\ < s, 
and w satisfies ||w|| L00 ( R d) < R, then 

\\x a (sdf(F(v + V )-F(v))\\ L2m <C £ \\x a (edf4 L2(Rd) . 

\a\+\0\<s M + l/3|<s 

Sketch of the proof of Lemma \4-b] When Xj? is replaced by (remove the control 
of the momenta), the result is exactly |^| Lemma 8.1]. The idea is to factor out w 
in the quantity F(y + w) — F(v) using the fundamental theorem of calculus, then to 
use Leibniz' rule, to conclude with Gagliardo-Nirenberg inequalities. In the case of 
Xg, the control of the momenta follows easily. □ 

We first notice that is uniformly bounded in L°°([Q, tq] x R d ). To prove that 
w £ N is bounded in L°°([0, To] X R d ), we use a continuity argument, and prove that it 
is actually small in that space, for N sufficiently large. This will be a consequence 
of the Gagliardo-Nirenberg inequalities: 

(4.5) for S >d/2, IMU-^) < \H H -m <e~ d/2 \Mx r 

(The scaling factor e~ d / 2 is obvious when one uses Fourier transform.) 
By construction, w|^(0, x) = O (e N+1 ) in any space X*. We first prove the result for 
N sufficiently large, then show how to get rid of this assumption. From Lemma l4~3*l 
there exists t(e, R) > such that if N + 1 > d/2, then for s sufficiently small, 

(4.6) \\v s N(t)h°°m < R 

for t G [0,t(e, R)]. As long as l|4.6|l holds, (|4.4I) and Lemma Wffl with s = imply 

dt h e N (t)\\ L * < C h e N {t)\\ L2 + Ce N \\r%(t)\\ L2 , 
and from Gronwall lemma, as long as (|4.ti|) holds for t < tq, we get that 

(4.7) \\^ N (t)\\ L2 <Cs N . 

The idea is now to obtain similar estimates for the momenta and derivatives of w^,. 
Applying the operator eV x to (|4.3|l yields: 

ied t (eV x v E N ) =H s {eV x v e N )+e\(t)(eV x ) (F(^) - F{v e N )) 
+ [eV,H s ] v e N -e N+1 eV x r e N . 
The same energy estimate as before gives: 

d t ||eV a w^(t)|| £a < HeV, (F(r) ~ F(v%))\\ L2 + \ ||[eV,H e ]^|| ia 



+ e N \\eV x r 



x' N\\ L 2 
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But we have 

[eV.H*] = (Wr)(j) + eVU{x). 
Since VVp is bounded and VU is sub-linear, the above estimate yields 

d t \\e\7 x v%(t)\\ L2 < \\eV x (F(^) - F(v e N ))\\ L2 + \ ||w^|| i2 + \\x^ N \\ L2 

(4 - 8) +s N \\eV x r%\\ L2 

< \\ey x ^ N \\ L2 + \\x^ N \\ L2 +e N - 1 , 

where we have used Proposition ^. 81 Lemma W7\ with s = 1, and (|4.7|) . We see that 
when £/ is quadratic, we have to find a similar estimate for || a^w^ Hz, 2 - For that, 
multiply i|4.3[) by x: 

ied t {xu e N ) = H e (xv s N ) + e\(t)x (F(yj £ ) - F(v%)) + [x, H £ ] - e N+1 xr%. 

Since [x,H e ] = —e 2 V x , the energy estimate yields, as long as (14.611 holds: 

5 t \\xu%(t)\\ L2 < \\x (F(r) - F(v%))\\ L2 + HeV^II^ + e N \\sV x r s N \\ L2 

<\\xv s N (t)\\ L2 + \\eV x v%\\ L2 +s N . 

Putting (|4.8jl and (|4.9|) together, we have: 

dt (||eV xW y £a + ||xw^(t)|| L2 ) < ||eV x w^[[ i2 + \\x^ N {t)\\ La + e N - 1 , 

and a Gronwall lemma yields, as long as (|4.6|l holds: 

(4.io) humx^^- 1 - 

One can check by induction that for k > 0, so long as Ij4.6|l holds, 
(4.H) h £ N (t)\\x S 

We now take advantage of the Gagliardo-Nirenberg inequality (|4.5fl . For s > d/2 
and as long as l|4.6J) holds, we get 

IKWIIl-^) < e- d/2 h e N (t)\\ Xi Z e N - s - d/2 ■ 

Thus, if N — s — d/2 > 0, a continuity argument shows that 14. 6|) holds up to time 
To provided that e is sufficiently small. In particular, w^-, hence ip s \ is well defined 
up to time To for < e < s(tq). To complete the proof of Theorem 14. 51 we have to 
prove (g21). Fix s, N G N; let s x > s such that s x > d/2, and Nj > sx +N+ 1. We 
infer from (|4.11|l that 

sup \\v £ Nl (t)\\ xS1 <e N i- s > <e N+1 . 
a<t<T„ 

It is straightforward that since Ni > N, 

sup llv^W-v^^ll^ <e N +K 

0<t<T 

We deduce that l|4.2jl holds for any s,N eN. □ 

Remark 4.7. A slightly shorter argument is available in the case d < 3, for which 
we have H 2 (R d ) C L°°(R d ), to prove Theorem in the case s = 2 only. The 
idea is to get an X 2 -estimate and use ()4.5[) again. Following an idea due initially 
to T. Kato consider the time derivative of the error w^. One can prove that 
\\edt^ £ N (t)\\L2 = 0(e N ), as long as (|4.6(l holds. Plugging this into (|4.3(l . we have, 
from (|4.7(l and since Vr is bounded and U is sub-quadratic: 

||e 2 Aw^(i)|| i2 <e N +\\x 2 ^ N (t)\\ L2 . 
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The control of ||x 2 w| vr (i)||^2 is then similar to (|4.9|) : 

\\x 2 u e N(t)\\ L2 < £ N + ||xMr(i)|| ia + ||e 2 Aw^(t)|| i2 , 
and we can conclude as above. 

Now it is easy to deduce the estimate announced in Theorem ll.il when ipj is as in 
Lemma 1X71 The L 2 estimate is (|4.2ll with N = s = 0. We have an L°° estimate, 
mimicking the above proof: for s > d/2 and N — d/2 > 1, 1)4.2)1 and 1)4.5)1 yield 

Sup \W(t)-V%(t)\\ Lx(Rd) <E- d ' 2 Sup \W(t)-V%(t)\\ X3 <E N - d l 2 <E. 

0<t<T O 0<t<T " 

It is straightforward that 

sup ||vg(t) -v s N (t)\\ LOO(Rd) <£, hence sup \\ip £ (t) - v £ Q {t)\\ Loo(Rd) < e . 

0<t<T O 0<t<T o 

Finally, we remove the assumption that ipf is as in Lemma 13.71 

Proposition 4.8. Let ip £ be the solution to (|l.l(l with initial datum tpj- as in 
Lemma \ tj. 7) Let ipj satisfying Assumptions \2.!A \2.1\ and Iff. 51 with K > d, and 
let ip £ be the solution to (|l.lf) with initial datum ipj. Then for any tq g]0, t[, there 
exists Bo > such that for < e < £q, ip e * s defined up to time tq. Moreover, 



sup 

0<t<T O 



ip £ (t)-i> e (t) =0(e K+1 - s ), fors>0. 



x; 



Remark 4.9. We deduce that Theorem 11.11 holds with an O (e d ) corrector in the 
initial datum. The L 2 estimate in Theorem 11.11 is straighforward, using Theo- 
rem 231 The L°° estimate Ijl.l9|l follows the same way, from 1)4. 5[l . Notice that 
the larger K, the more precise asymptotics we infer; for example, if K > d, we can 
remove the restriction r\ > in H1.19|l , using the above estimates and 14.5)1 . When 
s > K + 1, the above estimate does not look so good, since from Theorem 14.51 ip £ 
is bounded in X*. Yet, it gives some non-obvious control on ijj £ . 

Sketch of the proof of Proposition \4-8j The proof is very similar to that of Theo- 
rem 231 so we shall be brief. Introduce w e — ip £ — ip £ . It solves 

ied t v £ = H £ v £ + e\(t) (\ip £ \ 2a ij £ - |^ £ | 2<T V' £ ) , 

w £ | i=Q = O (e K+1 ) in X s e for any s G N . 

We can then follow the same lines as in the proof of Theorem l4.5l there is no source 
term (r% has disappeared) , and the size of w e is determined by the size of its initial 
datum. We have 

ll4=ollL~(R d ) £ e~ d/2 |Kt=olU| ^ e K+1 - d/2 , provided that s > | . 

Since K + 1 > d/2, we can start the argument of Theorem 14.51 Theorem 14.51 
and Sobolcv inequalities provide all the estimates we need for the "approximate" 
solution ip £ ; resuming all the arguments yields, so long as 1)4.6)1 holds, 

\\* B V)\\x;<e K+1 -'. 

Note that even if K + 1 — s < 0, we can apply a Gronwall argument to prove the 
above estimate. Since K + 1 > d, we can choose s > d/2 (not necessarily an integer, 
but this causes no trouble, by interpolation) such that K + l — s > d/2. The above 
estimate and 1)4.5)1 show that 1)4.6)1 holds up to time Tq, for e<1. □ 
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5. Generalizations and consequences 

5.1. Eigenvalue with multiplicity. As a first generalization we remark that all 
given results could be extended to the case where E n (k) is an isolated but m-fold 
degenerate family of eigenvalues, i.e. 

E n {k) = E*{k), Vne/CN, \I\ =m. 

Under the assumption (see e.g. [33] for a discussion on this) that there exists a 
smooth orthonormal basis {x;(fc, y))}i^j of ranll/(fc), where 

m 

nz(fc) :=X>(*)>fo(*OI 

1=1 

denotes the spectral projector corresponding to E*(k), the appropriate two-scale 
WKB-ansatz would then be 

m 

(5.1) ^ (t, x, Z ) ~ a °> l & (~ ' V ^(*> *)) e#(t,a;)/£ + 

2 = 1 

with (/>(i, a;) given by the solution of the Hamilton-Jacobi equation l|2.9|) with 
E n (k) — E*(k). As in |34l I45 j this would then lead to matrix-valued transport 
equations, which in our case are all coupled through the nonlinear term. The anal- 
ysis of this system is analogous to the scalar case but leads to rather intricate and 
tedious computations, which is why we neglected this situation. Also, from the 
physical point of view it is known that for periodic potentials such degeneracies are 
rather exceptional. (For the study of a similar 2-fold degenerated situation we refer 
to [42| . where a semi-classical scaled nonlinear Dirac equation is analyzed.) 

5.2. Wigner measures. Since Theorem 14.51 yields strong asymptotics for the 
wave-function in L 2 (R d ), we can compute the Wigner measure associated to the 
family ("0 £ )o<e<i- The Wigner measure of a family (ip e (t, '))o<e<i bounded in 
L 2 (M. d ) is the weak limit (up to the extraction of a subsequence) of its Wigner 
transform, 

(5.2) W E (t)] = J r (t, x- e -r,)W (t, * + |tj) e*" 

This limit is then found to be a nonnegative Radon measure on phase space. The 
Wigner transform has proved to be an efficient tool in the study of semi-classical 
and homogenization limits (see e.g. [Tl 1191 l2TH \'M\ ). 

Corollary 5.1. Let ip E (t) be the unique local-time-solution of (|1.1|) on [0,tq], as 
guaranteed by Theorem \4-5[ and let W e [ip e (t)] be its Wigner transform. Then, up 
to extraction of subsequences, we have 

(5.3) lim W e [^ e ] = fi in «S'([0, r ) x x Rl) weak-*, 

e—>0 q 

where the Wigner measure /i(t) ofip E (t) is given by 

2 

*(f-*-7*), 



(5.4) ^t,x,$= ] -j^dx< >^ 



7 *er* 



with k = V x 4>(t, x) S B. 
Proof. We have to compute 



lim f(x,0W e [r(t)](x ) 0dx<%= / f(x,0Kt,dx,d?), 



R 2./ 
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for any smooth test-function (observable) / £ S(W£ x R|). To this end, we plug 
the approximation Vg into the left hand side of this relation (that is, we use the 
strong L 2 convergence stated in Theorem II .10 . Since Xn(y,k) is T-periodic w.r.t. 
y £ M. d , we can rewrite it in form of a Fourier series: 

Xn(y,k) = -±- i £ e^' [ X n(z,k)e' lz ^dz. 

Using this representation, a non-stationary phase argument shows that all "non- 
diagonal" terms in (|5.2(l vanish in the limit e — > and hence (|5.4|1 is obtained from 
a straightforward computation. □ 

In our case, the strong convergence stated in Theorem 14.51 shows that the Wigner 
measure of (i^ s (t, 0)o<e<i ^ s tne same as m the linear case (see [201 Sect. 5.1]), since 
the main nonlinear effect appears as an order 0(1) phase oj, defined in Theorem ll.il 
In other words, the Wigner measure does not "see" the nonlinearity. This can be 
compared with the Wigner measures studied in [7|, for equations similar to Ijl.ip . 
without potential. For the same scaling as in (fTTTft . the main nonlinear effect was a 
"slowly" varying phase, which was invisible to the Wigner measure. It only appears 
as the first order correction in the Wigner transform. 

5.3. Complex-valued coupling factor. When the coupling factor X(t) is not 
real-valued, the analysis may be completely different; the approximate solution 
may blow up before the caustic. The first hint is that the L 2 -norm of ip e is not 
formally conserved. Multiply l|l.lll by ip e , integrate over R d and take the imaginary 
part: 

d 
dl 

On the other hand, the formal analysis of Sections Eland 01 still yields the transport 
equation (|3.4f> . which can also be written as l|3.8fl . Multiply (|3.8|l by So and take 
the real part: 

-\a (X t )\ =-lm K (X t ) lJf]a 

|2<7+2 



\\r(t)f L2 =21mX(t)\\r(t)\\T^ 



Y 

The solution of this ordinary differential equation may blow up in finite time before 
a caustic is formed, and the WKB-analysis breaks down at blow-up time. The 
above equation for the evolution of HV^W Hi 2 suggests that the exact solution may 
also blow up. In that case, the limitation for the validity of the WKB-expansion 
would not be a drawback of the method (as it is in the case of caustics), but a 
genuine nonlinear effect. 

Appendix A. Derivation of the leading order transport equation 

For the benefit of the reader, we shall discuss here in more detail how to pass from 
l|PJl to fOll . 

First, it will be convenient to rewrite (|3.2|) in a more symmetric form 
LxUq = id t u - - [D x ■ (D y + V x (f>) + (D y + V x <j)) ■ D x ] u Q , 
where from now on D x := —iS/ x . Then, inserting 

uo(t, x, y) = a (t, x)xn(y, Vx4>), 
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and denoting 

g n (t, x, y) = Xn (y, V x (j)(t, x)) , 
the solvability condition (|3.3|) can be written as 

d t a + (g n , d t g n ) L 2 (Y ) 0-0 + 7: (g n , V x • (D y + \7 x( j)) (a g n )) L2(Y) 
(A.l) , ' 

+ 2 (dn, (D y + V x cf>) • V a; (aog„)} i 2 ( y) - in(t, x)|a | 2cr ao = 0. 

Here we have used definition l|3.7(l and the fact that (xn,Xn) l 2 (y) = Differenti- 
ating the eigenvalue equation l|2.10(l w.r.t. to k yields 

(A.2) (V k H r (k) - V k E n (k)) Xn + (H T (k) - E n (k))V kXn = 0. 

Taking in this identity the scalar product with Xn we obtain 

, , . (Xn, V k H r (k)Xn} L 2 (Y) = (Xn, (Dy + k) X n} L 2 (Y) 

(A. 3) , , y ' 

= V k E n {k), 

since Hp is self-adjoint. From (|A.3() we deduce that (|A.lll can be written as 

(A. 4) d t a + {g n , d t g n ) L 2( Y ) a + V k E n (V x (f)) ■\' x a + f(t, x) a a = in(t, x)|a | 2<T a , 

where 

f(t,x) = \ ^ 9 "' ( D V + ' ^xg n ) L 2 {Y) + ^ (g n , V K • (D y + Vx4>)9n) L t(y) ■ 

Next, we substitute Xn by g n in (|A.3(I and differentiate w.r.t. x € M. d : 
(V x 5 n , (D y + V x (/))g n ) L2 ( Y) + (g n , V x ■ (D y + V x (j))g n ) L2{Y) = Aiv x V k E n (V x (j>). 
Since D y is self-adjoint and V x ^> is real, we have 

a ■= (9n, (D y + V x 4>) ■ V^g„) L2(y) = ((D y + V x 4>)g n , V x g„} i2( . y) , 
and we infer from above that 

a + A x 4> + a = div x V k E n (V x 4>). 

Therefore 

(A. 5) f(t,x) = a + ^A x (j) = Rea + ^A x (j) + ilma = i div^ V k E n (V x 4>) + i Im a. 



We simplify the last term. From (|A. 2|) . with k = \7 x <p, we obtain 

((D y + V x 4>) - V k E n (V x 4>)) g n + (H r (V x (t>) - E n (V x ct>)) V kXn (y, ^x<t>) = 0. 
Taking the L 2 (F)-scalar product by 

d 

dx 3 g n = ^x,ftx„ (y, V x <t>) 
i=i 

and taking the imaginary part, we have, since (%„, V^Xn) l 2 (y) e 
Ima = - iV ' k E n (V ' x (j>) ■ (g n ^x9n) L 2(Y) 

iAA " -^Im((ff r (V^)-S n (V^)) d kjXn , Y, 9 "^ 9 ^)- 

■j=i i=i 

The last sum also reads: 

J2 d 2 . Xi 0Im((if r (V x 0) - B„(V x 0)) d kiXn , d klX n) ■ 

l<j,l<d 
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Since Hr is self- adjoint, this term is zero. Hence, l|A.4jl together with (|A.5jl and 
(|A.6J) give the following equation for the principal amplitude: 

d t ao + (g n , dtg n ) L 2( Y) a Q + Ca Q + V k E n (V x (f>) ■ (g n , V x g n ) a Q = in(t, x)|a | 2<T a , 

where C is defined as in (|3.5|l . Finally, using the Hamilton- Jacobi equation H2.9J) , a 
straightforward calculation shows 

(A. 7) (g n , d t gn) L 2 (Y ) + ^kE n (V x ^) • (g n , V x g n ) = -f3(t,x) 

and we conclude that ao satisfies the nonlinear transport equation (|3.4|> . 

Acknowledgments. The authors are grateful to Stefan Teufel for fruitful discus- 
sions on this work. 
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